/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType hsTwophaseMixtureThermo<reactingMixture<gasThermoPhysics>>;

chemistryReader foamChemistryReader;

foamChemistryFile "$FOAM_CASE/constant/reactions";

foamChemistryThermoFile "$FOAM_CASE/constant/thermo";

// Phase volume fraction tolerance for diffusion mask
phaseMaskTol 0.01;

// Amount to clip alpha below to supress floatsom and jetsom.
phaseClipTol 1e-4;

Vapor
{
   Sc  Sc [ 0 0 0 0 0 0 0] 0.9;
   
   subspecies
   (
        CH2          {}
        CH2(S)       {}
        CH2O         {}
        CH2OH        {}
        CH3          {}
        CH3NHNH2     {}
        CH3NN        {}
        CH3NNH       {}
        CH3NNH2      {}
        CH3O         {}
        CH3OH        {}
        CO           {}
        CO2          {}
        H            {}
        H2           {}
        H2O          {}
        HCO          {}
        HNCO         {}
        HNO          {}
        HNO3         {}
        HONO         {}
        N            {}
        N2           {}
        N2O4         {}
        NCO          {}
        NO           {}
        NO2          {}
        O            {}
        OH           {}
   );
}

noVaporPairs
(
);


Liquid
{
   Sc  Sc [ 0 0 0 0 0 0 0] 200.0;
   
   subspecies
   (
        CH3NHNH2L
        {
            transportModel  Newtonian;
            nu              nu [ 0 2 -1 0 0 0 0 ] 9.2e-7; //rocketdyne
   
            D0           D0   [ 0 2 -1 0 0 0 0 ] 1e-9;
            rho0         rho0 [ 1 -3 0 0 0 0 0 ] 870.0;
            Tc           Tc   [ 0 0 0 1 0 0 0] 585.0;
            kappa        kappa [1 1 -3 -1 0 0 0] 0.25; //mmh bible, 300 K 
            
            sigma0          sigma0 [ 1 0 -2 0 0 0 0] 0.0626;
            sigmaa          0.884;
        }
   );
}



phaseChangeReactions
(    
    MMHEvap
    {
        phaseChangeModel LangmuirEvaporation;
        
        liquidSpecies
        (
            CH3NHNH2L   1.0
        );
        
        gasSpecies
        (
            CH3NHNH2   1.0
        );
        
        LangmuirEvaporationCoeffs
        {
            Lb              Lb [ 1 2 -2 0 -1 0 0] 38884000.0;
            Tb              Tb [ 0 0 0 1 0 0 0] 360.476;
            La              0.38;
            betaV           0.6;
            betaC           0.0;

            Tc              Tc [ 0 0 0 1 0 0 0] 585.0;
            Pc              Pc [ 1 -1 -2 0 0 0 0] 8210000.0;
            PvCoeffs        (9.1834 1943.89 492587.8);
        }
    }
);




// ************************************************************************* //
